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Verfahren zur Stabilisierung von technischen Prozessen 

Die Erfindung bezieht sich auf ein Verfahren zur Stabilisierung von technischen Prozes- 
sen resp. technischen Aniagen und ist geeignet zur Stabilisierung eiher verfahrens- und 
fertigungstechnischen Anlage. 

Es ist allgemein bekannt, fur die Qualitatsgrossen eines Prozesses Toleranzbander vor- 
zugeben, in denen ein Produktionsprozess ein gewunschtes Produkt produzieren soil. 

Weiterhin ist bekannt, wichtige Qualitatsgrossen eines Prozesses durch Anwendung 
mathematischer, statistischer oder neuronaler Algorithmen zu modellieren und vorher- 
zusagen, urn damit im voraus die Auswirkungen von Stellgrossenveranderungen auf 
den Produktionsprozess zu simulieren und nur bei gewunschten Resultaten fur die Qua- 
litatsgrossen die Einstellungen am Prozess real durchzufuhren. 

Es ist ebenso bekannt, mittels spezieller mathematischer Optimierungsverfahren, die 
optimalen Werte der wichtigsten Einflussgrossen fur einen vorgegebenen Qualitatsziel- 
wertes zu berechnen, siehe E, Schoneburg et al., Genetische Algorithmen und Evoiuti- 
onsstrategien, Addison-Wesley, Deutschland, 1994, S. 286-291 und S. 366-378 und US 
6.314.413 B1. 

Dennoch gibt es in technischen Prozessen sehr haufig Falle, wo der vorgegebene nu- 
merische Qualitatswert der Zielgrosse nur mit sehr grossen Aufwand oder auch gar 
nicht stabil realisiert werden kann, da der Produktionsprozess und mit ihm die Einfluss- 
grossen stets gewissen Schwankungen bzw. Storungen unterliegen. Deshalb gibt es 
sehr viele Produktionsprozesse bei dem das zu produzierende Produkt nicht mit ausrei- 
chender Stabilitat in der vorgegebenen Qualitat hergestellt werden kann. 
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Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren zur Stabilisierung von techni- 
schen Prozessen anzugeben, das fur die Qualitatsmerkmale eines technischen Produk- 
tes innerhalb der vorgegebenen Toleranzgrenzen denjenigen numerischen Zielwert 
vorgeben kann, bei dem dieses Produkt maximal stabil (oder auch flexibel) produziert 
wird und fur den zusatzlich alle einzusteilenden Werte der jeweiligen Steligrossen vor- 
geben werden konnen. 

Diese Aufgabe wird gelost durch ein im Anspruch 1 angegebenes Verfaiiren zur Opti- 
mierung der Produktion, bei dem mittels einer Global Scaling Kettenbruclianalyse die 
superflexiblen und superstabilen Werte der Qualitats- und Einflussgrossen berechnet 
und dem Produktionsprozess vorgegeben werden konnen. 

Vorteilhafte Ausgestaltungen sind in weiteren Anspruchen angegeben. 

Global Scaling (GS) ist ein eingefuhrter physikalischer Begriff, der verdeutlicht, dass die 
Verhaltnisse physikalisclier Grossen wie Massen, Temperaturen. Gewichte bei realen 
Systemen grosseninvariant sind und sich regelmassig logarithmisch wiederholen, siehe 
H. Muller, Global Scaling, Special 1, Ehlers Verlag 2001, S. 161/162. 

Erstmais aufgezeigt wurde dieser grosseninvariante Zusammenhang in biologischen 
Prozessen von Cislenko, der die Verteilung der Artenvielfalt der Flora und Fauna in Ab- 
hangigkeit ihrer Grosse und des Gewichtes uber eine logarithmisclie Gerade darstellte, 
siehe L. Cislenko, Die Stmktur der Fauna und Flora im Zusammenhang mit der Korper- 
grosse der Organismen, Moskau, 1981, S. 89-98. Dabei wurde gezeigt, dass die biolo- 
gischen Systems nicht beliebig auf der logarithmischen Grossenachse verteiit sind, 
sondern ihr Vorkommen in stets gleichen periodischen Abstanden auf der logarithmi- 
schen Geraden ihr Maximum bzw. Minimum haben. In technischen Systemen aussert 
sich dieser Zusammenhang darin, dass in einem technischen Prozess niemals alle 
durch den Prozess moglichen Prozesswerte gleichwahrscheinlich eingenommen wer- 
den konnen. 

Betrachtet man das Histogramm einer beliebigen technischen Messgrosse eines Inter- 
valls Zu bis Zo ungeglattet, so erkennt man neben der prozessabhangigen Hauptstruktur 
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des aufgenommenen Histogrammes, z.B. einer Gauss-, Poisson- oder Gleichvertei- 
lung, dass in der Feinstruktur dieses Histogrammes gewisse Messwerte haufiger auftre- 
ten, andere weniger haufiger. Mit bisher bestehenden Verfahren werden diese Uneben-. 
heiten weggefiltert, da sie als zufallige Storungen interpretiert werden. In der GS- 
Theorie gelten diese verschiedenen Auspragungen einer Messgrosse als wesentlich, da 
sie reproduzierbar sind. 

Mit Hilfe des OS lassen sich damit insbesondere diejenigen physikalischen Werte be- 
rechnen. die durch technische Prozesse bevorzugt eingenommen werden, denn nach 
OS konnen die von einem technischen oder naturlichen Prozess bevorzugten Werte 
durch eine Kettenbruclizerlegung nach Leonard Euler, Ober Kettenbruche,1737, Leo- 
nard Euler, Ober Schwingungen einer Saite, 1748, L. Euler, Leonhard Euler und Chris- 
tian Goldbach: Briefwechsel 1729 - 1764 (eds.: A. P. Juskevic, E. Winter), Abh. Deut- 
sche Akad. Wiss. Berlin, Akademie-Verlag 1965, ermittelt werden. denn nach Euler ist 
bekannt, dass jede reelle Zahl x durch Ihren Kettenbruch entsprechend Gleichung (1) 
dargestellt werden kann: 

X = no + z/ (ni + z/ (n2 + z/ (na + z/ (n4 + z/ (ns + ..) )))) (1) 

Die Grosse z stellt dabei den sog. Teilzahler dar, dessen Wert nach GS fur technische 
Optimierungen auf den Wert 2 festgelegt wird. 

Da die Skaleninvarianz in logarithmischen Massstaben auftritt, siehe Cislenko, werden 
im GS-Verfahren alle Anaiysen von zur Basis e iogarithmierten Grossen durchgefOhrt. 
Damit entsteht Gleichung (2) 

In X = no + 2 / (m + 2 / (nz + 2 / (ns + 2 / (n4 + 2 / (ns + ..) )))) (2) 

Die jeweiligen Zahlenwerte hangen von den zugrundeliegenden Masseinheiten ab, 
deshalb sind in GS fur alle gangigen Masse sog. Eichmasse y eingefQhrt, zu denen die 
auszuwertenden Grossen ins Verhaltnis gesetzt werden mussen. Dadurch entsteht die 
Gleichung (3) als spezielle Grundgleichung des GS, siehe H. MQIIer, Global Scaling, 
Special 1. Ehlers Veriag 2001, S. 157: 
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In (x/y) = no + 2 / (m + 2 / (02 + 2 / (n3 + 2 / (n4 + 2 / (05 + )))) (3) 

In Anwendungen von GS wird diese Grundgleichung (3) durch einen Winkel erweitert. 
mit (|) = 0 bzw. (j) = 3/2, urn den das logarithmische Verhaltnis In (x/y) vor eiher Ketten- 
bruchzerlegung verschobenen werden kann. 

Dadurch entstehtdie allgemeine Grundgleichung (4) des Global Scaling, nach der jede 
beliebige technische Messgrosse x zerlegt werden kanh: 

In (x/y) - (j) = no + 2 / (ni + 2 / (nz + 2 / (na + 2 / (n4 + 2 / (ns + ..) )))) (4) 

mit X ist die technische Messgrosse, gennessen in ihrer jeweiligen Einheit, y ist das na- 
turliche Eichmass dieser Grosse. Und <|)= 0 bzw. (|)= 3/2. 

Die Koeffizienten [no,ni,n2... ] mussen aufgrund der Konvergenzbedingung fur Ketten- 
bruche ihrem absoluten Betrag nach stets grfisser als der Zahler sein, siehe O. Perron. 
Die Lehre von den KettenbrQchen, TeubnerVerlag Leipzig, 1950, S. 62 und sind stets 
durch 3 teilbare ganze Zahlen. 

Diese Koeffizienten entscheiden Qber die charakteristischen Eigenschaften der Mess- 
grosse X, deshalb stellen sie einen sog. Kettenbruch-Code dar, Messgrossen deren 
Kettenbruchzerlegung nur einen Wert fur [hq] ergibt. befinden sich in einem Hauptkno- 
ten (Knoten .der Ebene no), Messgrossen mit Werteh von [hq, hi] in einem Subknoten 
der Ebene ni usw. Kernbereich eines Knotens n\ ist der Bereich in unmittelbarer Nahe 
eines Knotens. d.h. Subknoten nj liegt im Bereich -9 bis 9. mit j=i+1. Nach der GS- 
Theorie sind nun folgende physikalischen Eigenschaften in Abhangigkeit des Ketten- 
bruch-Codes [no. ni.n2, na...] bekannt: 
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Kettenbruch-Code Beispiel 


Eigenschaft 


Bemerkungen 


[no +3 -3 +3 -3,..] 


Superflexibilitat 


Messgrosse befindet sich am Rand eines 
Knotenbereiches, in maximaler Entfer- 
nung vom stark fluktulerenden Kernbe- 
reich, in elnem sog. superflexiblen Be- 
reich. 

Das zugrundeliegende System reagiert 
flexibel auf Verandeaingen und StSrun- 
gen 


[no 54...] 
[no 57...] 
[Ho 60...] 


SuperstabilitSt 


Messgrosse befindet sich auf einem Kno- 
ten. ist deshalb extrem stabil. Kann aber 
nur bei genauer Prazision erreicht wer- 
den. ansonsten ist MessgrSsse im stark 
fluktulerenden Kembereich des Knotens 


[noni... 4]... [noni... 5] 
[no Ri... 7].,. [no Hi... 8] 


Relative Ruhe 


Messgrosse befindet sich zwischen zwei 

Subknotenbereichen. 

Messgrosse ist damit in einem Berelch 

relaliver Fluktuationsarmut 



Damit kann fur jede beiiebige technische Messgrosse nach einer GS-Analyse erkannt 
werden, ob diese Grosse stabil ist. grossen Schwankungen unterliegt. flexibel auf Sto- 
rungen reagiert oder in einem Bereich relativer Fluktuationsarmut ist. Alle weiteren 
Kombinationen ergeben sich aus diesen Uberlegungen. 



Beim erfindungsgemassen Verfahren wird das Ziel verfoigt, nicht nur die optimalen 
Werte des Produktes, z.B. das Drehmoment eines Maschinenteils, die optimale Glanz- 
charakteristik eines Lackes oder die optimale Haftung eines Reifens zu ermitteln, son- 
dern gleichzeitig diejenigen Werte der Einflussgrossen zu berechnen, die letztendlich zu 
dem jeweiligen Produktwert fiihren und die gleichzeitig einen besonders robusten, stabi- 
len bzw. flexiblen Produktionsprozess bedingen. 

Die Aufgabe wird erfindungsgemass mit folgenden Verfahrensschritten gelost: 

1 . Bestimmung der vorhandenen Qualitatsgrossen und ihrer zulassigen Toleranzen 
(Eingabe der Toleranzen fur Zielgrossen) 

2. Global Scaling Analyse innerhalb der vorgegebenen Toleranzen fur jede Quall- 
tatsgrosse zum Ermitteln aller GS-optimalen Werte dieser Qualitatsgrosse 

3. Aufnahme der Prozessdaten 

4. Prozessmodellierung und Sensitivitatsanaiyse der vorhandenen Einflusse zur Er- 
mittlung der lokalen Prozesseinflusse, die letztendlich zur Produktion mit diesem 
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Qualitatswert gefuhrt haben und damit zur Bestimmung der wichtigsten Einfluss- 
grossen 

5. Prozessoptimierung fur den bzw. fur alle vorgegebenen GS-optima!en Zielgros- 
senwerte durch Ruckrechnung auf die Werte der wichtigsten Einflussgrossen 

6. Durchfuhrung einer Global Scaling Analyse der wichtigsten Einflussgrossen fur die 
Auswahl des GS optimalen Zielgrossenwertes, wenn es aus Produktsicht im vor- 
gegebenen Intervall mehrere gleichberechtigte GS-Werte fur das Produkt gibt. 

7. Bestimmung des optimalen Prozesszustandes und Ruckrechnung der dazugeho- 
rigen Einflussgrossen, insbesondere der Stellgrossen (Auswahl des GS-optimalen 
Zielgrossenwertes) 

8. Ausgabe der Stellgrossen an den Prozess und ggf. Ruckkoppiung zu Schritt 3 
(und RQckrechnung) 

Jeder einzelne Verfahrensschritt wird im weiteren detailliert am Beispiel der Herstellung 
von Scharnieren eines Automobilzulieferers an Hand einer Zeichnung beschrieben: 

in der Zeichnung zeigen die: 

Fig. 1 Gemessenes Drehmoment MD eines produzierten Scharnieres uber die Zeit 
(Zeitachse in 10 Sekunden-Takt) 

Fig. 2 Histogramm des Drehmomentes MD mit dem Messbereich O Nm bis 2,5 Nm 
auf der X-Achse 

Fig. 3 Bin Histogramm des Drehmomentes MD mit dem Messbereich 1,0 Nm bis 
1 ,5 NM auf der X-Achse 

Fig. 4 Feinstruktur von Histogrammen zweier Produktionslinien fur rechte und linke 
' Scharniere. die 3 Tage lang zeitgleich produziert wurden 

a), Scharniere links MDL b) Scharniere rechts MDR 



Fig. 5 



Werkzeug GSC3000 zur GS-Analyse von physikalischen und technischen 
GrOssen 
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GS-Analyse des Drehmomentes (In Ncm) von Kfz-Scharnieren 

Histogramm der erzeugten Drehnnonnente im Bereich 0,87 bis 0,92 Nm 

Sensitivatsanalyse einer Scharnierproduktion (Sensitivity der Inputs zur Out- 
put-Variablen MD) 

Prozessoptimierung mitels Self-Organizing Maps durch automatische Ruck- 
rechnung von der Zlelgrosse MD auf die einzustsellenden Eingangswerte 

Darstellung von mogllchen Eingangsgrossen fiir Scharniere mit Zielgrosse 
Drehmoment MD = 1 ,392 Nm . 

Histogramm derTemperatur TW einer Produktlonslinie fur die Fertigung von 
Schamieren 

1 . Aufnahme der Qualitatsgrossen und der zulassigen Toleranzen 
Die wesentliche Qualitatsgrosse eines Scharniers ist neben seinen geometrischen Ab- 
massen das sogenannte Drehmoment oder Schwenkmoment MD. Durch die Anwen- 
dungen des Scharniers im hier dargestellten Anwendungsfall muss das Drehmoment in 
einem engen Toleranzband, im dargestellten Beispiel 0,50 - 2 Nm, produziert werden. 
Welche Werte es innerhalb des Toleranzbands einnehmen soli, ist durch die Anwen- 
dung nicht explizit vorgegeben. Im vorliegenden Montageprozess ist durch eine 100%- 
Qualitatskontrolle nach der Montage fQr jedes einzelne Scharnier das Drehmoment MD 
ausgemessen worden und beispieihaft Qber einen Zeitraum von mehreren Tagen in Fig. 
1 dargestellt: 

Die Figuren 2 und 3 stellen das Histogramm des aufgenommenen Drehmomentes MD 
mit der Inten/allbreite 1 Ncm einmal fUr den gesamten Messbereich von 0 -2,50 Nm 
und einmal fQr den Teiibereich 1 ,0 - 1 ,50 Nm dar. 



Fig. 6 
Fig. 7 
Fig. 8 

Fig. 9 

Fig. 10 

Fig. 11 
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Die Feinstrukturder Histogramme in beiden Darstellungen zeigt die unterschiedliche 
Auspragung der realisierten Drehmomentwerte im realen Prozess. Beispielsweise wird 
derWert 1.26 Nm (126 Ncm) im dargestellten Zeitraum 70 Mai, der Wert 1,27 Nm (127 
Ncm) dagegen nur 30 Mai realisiert. Die nachfolgende Tabelle gibt elnen Oberblick u.ber . 
die Extremwerte dieser Haufigkeitsverteilung: 



■ GemesseiiaiiS 


; Extnemwerttyp- /.r;.";^^:^":: 


; Mittelgyh|<f| dej;jT^ 


>64 


Haufungslntervalle 


101, 113.126.130. 134 


<42 


LOckenintervaile 


100, 104. 106, 112. 127 



Diese Feinstruktur 1st kein Zufall, da sie sich bei jedem anderen Zeitraum in einer ahnli- 
chen Struktur stets reproduzieren lasst und fur zeitgleich produzierte Scharniere sogaf 
ahnlich ist wie Fig. 4 zeigt. 

Man erkennt, dass fur beide Produktionslinien eine Feinstrukturder Histogramme exis- 
tiert und diese eine gewisse Ahnlichkeit besitzen. Untersucht man die Produktionslinien 
uber mehrere Monate so stellt man fest, dass es auf jeder Produktionslinie und an je- 
dem Tag Drehmomente gibt, die besonders haufig produziert werden und Drehmomen- 
te, die besonders selten produziert werden. Da die GS-Theorie genau eine solche Ver- 
teilung und Histogrammstruktur von physikalischen und technischen Messgrossen vor- 
aussagt und sich die Produktionslinien von Scharnieren damit GS-konform verhalten, 
wird eine GS-Optimierung durchgeftihrt, um eine GS-optimaie Produktionslinie aufzu- 
bauen. 

Im weiteren wird beispielhaft die GS-Analyse des Drehmomentes beschrieben, um die 
stabilen bzw. flexiblen Drehmomentenwerte zu berechnen. 

2. Global Scaling Analyse der Zielgrosse Drehmoment 
2.1 . Kettenbruchzerlegung des Drehmomentes 

Gegeben ist ein Drehmoment MD in einem Toleranzbereich zwischen 50 Ncm und 200 
Ncm. Gesucht ist der nach GS-optimaie Wert innerhalb des Toleranzbereiches. Das 
Eichmass y fur das Drehmoment ist 1,503277E-10 Nm = Protonenmasse * c^. siehe 
Absatz Eichmasse in H. Muller. Global Scaling. Speciall, Ehlers Verlag 2001. S. 129. 
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Nach Gleichung (4) mit ^ = 3/2 ergibt sich eine Kettenbruchzerlegung und die Berech- 
nung der Koeffizienten no. ni, na usw. Der Wertebereich 50 Ncm - 200 Ncm entspricht 
dem Knotenbereich I24±1]. also 50 Ncm = [24; -3] und 200 Ncm = [24; +3]. 

Die Berechnung der Drehmomentenwerte durch KettenbrQche nach Gleichung (4) wur- 
de beispielhaft .mit dem Werlczeug GSC3000 des Institutes fur Raum-Energie- 
Forschung; MQnchen (IREF) durchgefOhrt und ist in Fig. 5 exemplarisch fur den Dreh- 
momentwert MD = 1,27 Nm dargestellt. 

Die lineare Mitte des Toleranzbereiches 0.5 Nm und 200 Nm und vorgegebene Ziel- 
grosse einer Steuerung fur MD ist 125 Ncm = [24; +6. -3. +3, ...] was die Dominanz des 
Wertes 125 Ncm erklart. Der Wert 127 Ncm ist nach GS gleich [24; +6. -5]. In den . 
Randbereichen [24; -3] bzw. [24; +3] bzw. in den Sublucken. z.B. zwischen [24; +5] = 
133 Ncm und [24; +4] = 146 Ncm, sind die Fluktuationen minimal ausgepragt. Maximale 
Fluktuationen sind innerhalb des Kernbereichs von [24; -10] = 73 Ncm bis [24; +10] = 
109 Ncmzu enA^arten. 

Die graphische Darstellung in Fig. 6 der Drehmomentwerte in AbhSngigkeit von no und 
ni und der Phase ^ verdeutlichit die Zusammenhange fQr die Werte <|> = 0 bzw. ^ = 3/2, 
no=24 und ni =-3 bis -» bzw. ni =+3 bis +». Die durch den Kettenbruch nach Glei- 
chung (4) berechneten Werte des Drehmomentes von 33 Ncm bis 1082 Ncm sind dabei 
Qberdie LSngsachse der Fig. 6 aufgetragen. 

Ware der Prozess der Scharnierberstellung weder geregelt noch gesteuert, so wGrde 
sich die Haufigkeitdertechnisch produzierten Drehmomentwerte nach Fig. 6 verteilen. 
Die weissen Flachen in Fig. 6 verdeutiichen Drehmomentwerte, die relativ selten auftre- 
ten (sog. Luckenbereiciie), die grauen Fiachen verdeutiichen Bereiche mit hSufigen 
Drehmomenten und die dunkelgrauen Flachen stellen sog. Knotenbereiche des Dreh- 
momentes dar. 

Da aber im vorliegenden Produktionsfall der Produktionssteuerung der Zielwert 
IVtDsoLL=125 Ncm vorgegeben wurde. weichtdie Verteilung dertechnisch realisierten 
Drehmomente von dieser theoretischen Verteilung ab und stellt sich nach Fig. 2 dar. 
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Nach GS konnen Werte in den Knotenbereichen, hier um 89 Ncm, 124 Ncm und 173 
Ncm, technisch nicht stabil produziert werden, da in Knotennahe stets Fluktuationen 
und Storungen auftreten werden, siehe Tabelle. Seite 5 und exemplarisch Fig. 7. 

In Fig. 7 erkennt man, dass es um den Bereich 0,88 Nm bis 0,92 Nm zu grossen Fluk- 
tuationen dertatsachlich produzierten Drehmomente von Schamieren kommt, was be- 
deutet. dass der Wert 0.89 Nm in der Praxis nicht stabil produziert werden kann. 

FQr die Global Scaling Optimierung der Produktion von Schamieren gibt es im wesentli- 
chen zwei Ziele: 

Das erste Ziel besteht in der Einhaltung der zu produzierenden Produktqualitat. Das 
Drehmoment MD jedes Scharniers soil moglichst in der Mitte des Toleranzintervalls mit 
minimalen Toleranzschwankungen AMD produziert werden. Optimal sind deshalb 
Nenngrossenbereiche mit minimalen Fluktuationen, also Werte zwischen 

MDsubi = 133 Ncm = [24; +5] und MDsub2 = 146 Ncm = [24; +4], 

Der optimale Nennwert fur die Produktion von Schamieren im Toleranzband 50 bis 200 
Ncm liegtdamitin der logarithmischen Mitte von 133 Ncm und 146 Ncm. also be! 139,2 
Ncm. 

Das zweite Ziel be! der Produktion von Schamieren ist es, den Prozess der Scharnier- 
herstellung moglichst robust (also flexibel) gegenuber StOrungen zu realisieren. Das 
Drehmoment MD entsteht im Produktionsprozess im wesentlichen durch das sog. 
Obermass. Das Obermass ist der Bolzendurchmesser des Scharniers plus zweimal die 
Buchsenwandstarke. es ist grosser als Scharnierlochdurchmesser, in den der Bolzen 
gepresst wird. 

Die Einzelteile des Scharniers unterliegen in diesem Anwendungsbeispiel Toleranzen, 
die aus KostengrQnden, wegen der mechanischen Fertigung, nicht welter reduziert wer- 
den kSnnen. deshalb kommt es darauf an. die Nenngrossen auch der Einzelteile inner- 
halb jedes Toleranzbandes nach GS optimal zu bestimmen. Schamierbolzen und 
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Scharnierlochdurchmesser liegen im Bereich urn 12 mm, die Buchsenwandstarke bei 
0,48 mm. Urn die Abhangigkeiten des Drehmomentes MD von den Einzelteilen zu mo- 
dellieren, sind verschiedene empirische Modelliemngsverfahren bekannt. 

In diesem Anwendungsfall wird ein empirisches Data Mining Verfahren zur Modellierung 
beschrieben. da as auf grOsstmSgliche Prazision des funktionalen Zusammenhanges 
zwischen Drehmoment und EinflussgrCsseri ankommt. Jede andere Modellierung. zum 
Beispiel sine lineare Modellierung zwischen Obermass und Drehmoment funktioniert 
entsprechend. hat aber geringere Modellgenauigkeiten. 

3. Aufnahme der Prozessdaten 

Es ist Stand der Technik, dass mittels verschiedener Messaufnehmer und Leittechnik 
die Prozessdaten im vorgegebenen Prozesstakt, im Anwendungsbeispiel 10 Sekunden, 
mitgeschrieben und archiviert wierden. 

4. Prozessmodellierung und Sensitivitatsanalyse fOr die ZielgrSsse 

Es ist Stand der Technik zur Prozessmodellierung verschiedene Verfahren. analytische 
(wie Differentialgleichungen) oder empirische (wie lineare und nicht lineare Regressio- 
nen oder neuronale Netze einzusetzen. je nachdem wie gut der Prozess analytisch be- 
schrieben werden kann. FQrfertigungstechnische Prozesse. wie z.B. komplexe Monta- 
gen, werden auch wissensbasierte Verfahren verwendet, bei denen das Verhalten des 
Prozesses anhand seiner Prozessdaten gelernt wird, siehe M. Poike, Prozessleittech- 
nik, Oldenbou'rg verlag, Munchen/Wien 1994, S. 813-817. 

Im voriiegenden Beispiel wurde die Prozessmodellierung durch eine neuronale Netz- 
modellierung realisiert. Die erreichte Genauigkeit liegt bei knapp 88%, d.h. mit knapp 
88%iger Genauigkeit bezogen auf den Messbereich der Drehmomentmessung kann mit 
dem-Modell aus 10 EingangsgrSssen der Drehmomentwert geschatzt werden: 
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Es ist weiterhin bekannt, die quantitative Abhangigkeit einer Qualitatsgrosse, im vorlie- 
genden Beispiel das Drehmoment MD, von seinen Einflussgrossen durch eine modell- 
baslerte Sensitivitatsanalyse zu realisieren, siehe R. Otte. Selbstorganisierende Merk- 
maiskarten zur multivariaten Datenanalyse komplexer technischer Prozesse, Shaker 
Verlag, Aachen, 1999 und R. Otte et al.. Data Mining fur die industrielle Praxis, Carl 
Manser Verlag, 2004. Die Fig. 8 zeigt das Ergebnis einer Sensitivitatsanalyse basierend 
auf den oben genannten neuronalen Modell der Scharnierproduktion. 

Die Fig. 8 zeigt. dass die wichtigste Einflussgrosse fur das Drehmoment MD der Bolzen 
DB1. dann der Bolzen DB2. dann der Kalibrierdorn DK ist. Die Temperaturen TW und 
TU haben den geringsten Prozesseinfluss auf das Drehmoment. Das bedeutet, dass die 
Anderungen der Temperatur im Mittel zu relativ geringeren Andemngen des Drehmo- 
mentes fuhren werden. 

Durch den Verfahrenschritt 4 werden automatisch diejenigen Eingangsgrossen identifj- 
ziert, die einen Einfluss auf den Wert der Zielgrosse, im Beispiel das Drehmoment MD, 
haben. 

5. Prozessoptimierung fQr einen GS-optimalen Wert der Zielgr5sse 
Im Verfahrenschritt 5 werden fur die wichtigsten EinflussgrSssen bzw. Steilgrossen des 
Verfahrensschrittes 4 die konkreten quantitativen Nennwerte bestimmt, um den vorge- 
gebenen Zieigrossenwert zu erreichen. 

Es ist Stand der Technik empirisch gelernte Modelle fQr Aufgaben der Optimierung zu 
invertieren. um fOr gegebene Zielwerte, die dazugehorigen Werte der Einflussgrossen 
zu berechnen. Beispielsweise kann das Verhalten des Prozesses durch Self-Organizing 
Maps abgebildet werden, siehe EP 0 845 720 B1 und R. Otte, Selbstorganisierende 



wo 2005/071504 



PCT/CH2005/000013 



13 

Merkmalskarten zur multivariaten Datenanalyse komplexer technischer Prozesse, Sha- 
ker Verlag, Aachen, 1999. Die Fig. 9 zeigt eine derartige modellbasierte Datenanalyse, 
die auf dem Verfahren der Self-Organizing Maps beruht 

Durch Verwendung eines Prozessmodells zu Optimierungsaufgaben kann fur jeden 
vorgegebenen, GS-optimalen Wert der Zielgrosse, im Beispiel MD, berechnet werden, 
welche Eingangsgrossen letztendlich zu dem Zielwert fuhren. Dannit werden erfin- 
dungsgennass fur jede GS-Berechnung des Zielgrossenwertes die dazugehorigen Wer- 
te der Einfiussgrossen zuruckberechnet 

Die Fig. 10 zeigt beispielhaft die Werte ausgewahlter Eingangsgrossen Kalibrierdorn 
DK, Temperaturen, TW und TU, Scharnierhalftendurchnnesser DS, Buchsenwandstarke 
BW und Bolzendurchmesser DB1, DB2, DB3 fur die GS-optimale Zielgrosse MD = 
139.2 Ncm = 1,392 Nm. 

Damit sind die Eingangsgrossen bestimnrit, die zu einem GS-optimalen Zielgrossenwert 
Ygs, im Beispiel MD = 1,392 Ncm, fuhren. 

6. Global Scaling Analyse der wichtigsten Einfiussgrossen 

Durch Anwendung des Verfahrensschrittes 1 und 2 ist die Bestimmung aller GS optima- 
Ian Werte Ygsi der Zielgrosse im vorgegebenen Toleranzband ymin und ymaxdurchge- 
ftihrt. Durch Anwendung der Verfahrensschritte 3, 4 und 5 sind die dazugehorigen Ein- 
gangsgrossenwerte des Prozesses ermittelt. 

Wenn es im Intervall aus Produktsicht nur einen einzigen optimalen GS-Wert yes gibt, 
dann wird der hier beschriebene Verfahrensschritt 6 Qbersprungen und die zu dem GS- 
Wert yes zugehorigen StellgrSssen nach Verfahrensschritt 5 und 7 berechnet. Sind 
mehrere gleichberechtigte GS optimale Werte yosi vorhanden. dann folgt die Auswahl 
durch die GS-Analyse der Eingangsgrossen durch mehrmalige Anwendung des Verfah- 
rensschrittes 5. d.h. der RQckrechnung der zu einem Optimum yosi fuhrenden Ein- 
gangsgrossen. 
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Damit sind die Zahlenwerte der Eingangsgrossen bestimmt. Auch diese Zahlenwerte 
unterliegen der GS-Struktur, wie die Fig. 1 1 am Beispiel der Temperatur TW zeigt. 

Nach der Ermittlung der konkreten Eingangsgrossenwerte nach Verfahrensschritt 5 wird 
far jeden der gefundenen Eingangswerte Xni die GS-Analyse nach Verfahrensschritt 2 
durchgefuhrt. Damit werden fur die wichtigsten Eingangsgrossen die OS optimalen 
Werte ermittelt. Die GS- optimalen Eingangswerte hangen von dem Typ der Eingangs- 
grosse ab. 

Die Eingangsgrossen sind fur folgende Typen GS-optimal: 

a) Eingangsmessgrossen: Optimum nach GS durch Ermittlung der maximalen Auf- 
trittswahrscheinlichkeit. da diese Wertrealisierungen auch in Zukunftam haufigs- 
ten auftreten werden, 

b) Sind die Eingangsmessgrossen IVlesswerte von technischen Produkten so konnen 
die Nenngrossen fur diese Werte nach GS optimiert werden. 

Im vorliegenden Fall ist beispielsweise die Buchsenwandstarke BW im Bereich 
0,46 mm. bis 0,5 mm zu optimieren. 
BWunten = 0,46mm = [21 3 3 -3 3 5] 
BWNenn = 0,48mm = [21 3 3 21 6 -9] 
BWoben = 0,5mm = [21 3 6 --3 3 -3] 

Als Nenngrosse ist damit der Wert BW= 0.471 mm nach GS gleich [21 3 3 -5] und 
damit optimal, da er in einer Lucke auf na liegt. 

c) Eingangsstellgrossen: Optimierung nach GS hinsichtlich Superflexibilitat, da diese 
Werte am unempfindlichsten auf Anderungen im Wert selber und im Prozess 
reagieren. 

Im Verfahrenschritt 6 wurde fur jeden GS-optimalen ZielgrSssenwert. nach Verfahrens- 
schritt 2. und seiner dazugehorigen Eingangsgrossen. nach Verfahrensschritt 4 und 5, 
eine GS-Analyse durchgefuhrt, wodurch eine vollstandige GS-Analyse des Prozesses 
realisiert ist. 

7. Bestimmung des optimalen Prozesszustandes 
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Im Verfahrensschritt 7 wird aus der Menge aller nach GS moglichen ProzesszustSnden, 
der GS-optimale Prozesszustand ausgewahlt und die dazugehorigen Stellgrossen er- 
mittelt. 

Das Prozessoptimum ist fur die maximale Anzahl von GS-optimalen Eingangsgrossen 
erreicht, fur die die Zielgrosse ycsi gerade GS-optimal ist. Sind mehrere Zielgrossen- 
werte gleichberechtigt, so ist derjenige Wert ycsi der weiter in der logarithmischen Mitte 
des Intervalls ymin und Vmax liegt der optimale. 

Aus der Menge der zunn gefundenen Optinnum Ygs gehorenden Eingangsgrossen Xigs. 
X2GS. X3GS XnGs'wird die Untermenge der Stellgrossen xsteiiiGs ausgelesen. 

Damit ist das Optinnunn fur den fertigungs- oder verfahrenstechnischen Prozess gefun- ' 
den und die Stellgrossen xsteincs zur Erreichung des Optinriums werden tiber die vor- 
handene Leittechnik an den Prozess ausgegeben. 

8. Ausgabe der GS-optimalen Werte an den Prozess 

Die GS optimierten Werte der Stellgrossen des optimalen Prozesszustandes nach. 
Verfahrensschritt 7 werden beispielsweise uber eine Leittechnik oder eine SPS an den 
technischen oder chemischen Prozess ausgegeben, womit der Prozess 
erfindungsgemass nach GS optinniert ist. 
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Patentanspruche 

1 . Verfahren zur Stabilisierung von technischen Prozessen, gekennzeichnet durch 
die nachstehenden Schritte: 

a) Aufnahme und Auswahl der Daten und Aufbereitung der Daten, 

b) GS Analyse des Zielwertes innerhalb des vorgegebenen Toleranzbandes 

c) Automatische Prozessmodellierung und Sensitivitatsanalyse zur Bestimmung 
der fur die gewahlten Zielgrossen Y\ wichtigsten Einflussgrossen Xi, 

d) Auswahl der im Toleranzband GS optimalen ZielgrOssenwerte und Ruck- 
rechnung auf die dazugehorigen Eingangsgrossenwerte, Nennwerte, aller 
Einflussgrossen durch gangige Optimienjngsverfahren 

e) GS Analyse der unter d) berechneten Nennwerte der Einflussgrossen und 

f) Wiederhoiung der Verfahrensschritte d) und e) bis zu einem GS optimalen 
Zielgrossenwert aus dem vorgegebenen Toleranzband alle wichtigen GS op- 
timalen Nennwerte der Einflussgrossen und Stellgrossen gefunden wurden 

g) Ausgabe der GS optimalen Stellgrossen an den Prozess 

2. Verfahren zur Stabilisierung von technischen Prozessen, dadurch gekennzeichnet, 
dass mehrere Zielgrossen gleichzeitig und zusammenhangend nach Global Sea- 
ting (GS) optimiert werden, indem die ausgewahlten Zielgrossen durch eine belie- 
bige mathematische Funktion miteinanderverknQpft werden und die resultierende 
Gutefunktion (G) GS-optimiert wird. 

■ 3. Verfahren nach Anspmch 2, bei dem alle fur die Optimierung relevanten Grossen 
durch eine Analyse auf Grundlage von Global Scaling Verfahren hinsichtlich Fiexi- 
bilitat und Stabilitat optimiert werden, indem eine Global Scaling Analyse von Da- 
ten der relevanten Parameter und die dazugehorigen Zielgrossen realisiert wird. 
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4. Verfahren nach Anspruch 3, dadurch gekennzeichnet, dass der Prozess selbst 
optimiert wird. 

5. Verfahren nach Anspruch 3, dadurch gekennzeichnet, dass die Produkte die durch 
einen Prozess hergestellt werden optimiert werden. 

6. Verfahren nach einem der vorstehenden Anspruche 2 bis 5, dadurch gekenn- 
zeichnet, dass zur GS Optimierung einer Zielgrosse auch die Einflussgrossen des 
Prozesses einer GS Analyse und Optimierung unterzogen werden. 

7. Verfahren nach einem der Anspruche 2 bis 6, dadurch gekennzeichnet, dass Pro- 
zesse resp. Produkte des IVlaschinenbaus oder der chemischen Industrie stabili- 
siert bzw. optimiert werden. 
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Drehmoment MD uber die Zeit 




Zeitachse im 10 sec. Takt 

Fig. 1 Gemessenes Drehmoment WID eines produzierten Scharnieres uber die Zeit (Zeitachse 
in 10 Sekunden-Takt) 




Fig. 2 Histogramm des Drehmomentes MD mit dem Messbereich O Nm bis 2,5 Nm auf der X- 
Achse 
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Fig. 3 Ein Histogramm des Drehmomentes MD mit dem Messbereich 1,0 Nm bis 1,5 NM auf 
der X-Achse 
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Fig, 4 



Feinstruktur von Histogrammen zweier Produktionslinien fur rechte und linke Schar- 

niere, die 3 Tage lang zeitgleich produziert wurden 

a), Scharniere links WDL b) Scharniere rechts MDR 
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Fig. 5 Werkzeug GSC3000 zur GS-Analyse von physikalischen und technischen Grossen 



Global Scaling der Drehmomente (Ncm) 




^ -6-9 oo +9 +6 +3 



Fig. 6 



GS-Anaiyse des Drehmomentes (In Ncm) von Kfz-Scharnieren 
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Fig. 7 Histogramm der erzeugten Drehmomente im Bereich 0,87 bis 0,92 Nm 




Fig. 8 Sensitivatsanalyse einer Scharnlerproduktion (Sensitivity der Inputs zur Output 
Variablen MD) 
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Fig. 9 Prozessoptimierung mitels Self-Organizing Maps durch automatische RQckrechnung 
von der Zielgrosse MD auf die einzustsellenden Eingangswerte 



1234 ! AAB fi 1 0.019' 



2«ISB : RTCfR 



UKKIR 



' 0,019; 

" o,oi9r 



0.01 7 j 



90 



90 



80 



70 



90 



23.611 25.231 0.01382| D.OODdBSi 0,013957 




"^^ra! 2417] 0^401 1 DjgDjTH 0J013962 ! 0.013958! O .QI.^ 

■ T'- : --L^. -X^ ^.^^i 7r'A r^r,,ni, n flT^ci nni-acCHl n m-^QA-y 



' Z2.39I 



23.S5i 
24 71 



aQT^rOfflX )JB4| 6.013961 0.0135^ 1 0.013 942 ! 1.357 



"OJUOSi 00304681 o toi 3969! 0 pi 3S66i 0.0139431 0.04Gi 



a^ai 23^51 o5l313l ap004h; 0.013967'!" 0.01 396B! 6,013947} 1.37 



r435i 1 .392 



1.438 



L445 



1.419 



■ 1.392 



1.392 



1.392 



12 



1357 



AAH R I O^ i 



gP] 24361 7A2\\ 0 !o 1199 1 0J000471; q!o1 3 98 7 y Ofll 3867 ! 0.01 3955 j 



Fig. 10 Darstellung von moglichen Eingangsgrossen fiir Scharniere mit Zielgrosse Drehmo- 
mentMD = 1,392 Nm 




Fig. 11 



Histogramm der Temperatur TW einer Produktionslmie fiir die Fertigung von Schar- 
nieren 



